Evolutionary conservation of an ancient retroviral gagpol gene in Artiodactyla.

The genomes of mammals contain fingerprints of past infections by ancient retroviruses that invaded the germline of their ancestors. Most of these endogenous retroviruses (ERVs) contain only remnants of the original retrovirus; however, on rare occasions, ERV genes can be co-opted for a beneficial host function. While most studies of co-opted ERVs have focused on envelope genes, including the syncytins that function in placentation, there are examples of co-opted gag genes including one we recently discovered in simian primates. Here, we searched for other intact gag genes in non-primate mammalian lineages. We began by examining the genomes of extant camel species, which represent a basal lineage in the order Artiodactyla. This identified a gagpol gene with a large open reading frame (ORF) (>3,500 bp) in the same orthologous location in Artiodactyla species but that is absent in other mammals. Thus, this ERV was fixed in the common ancestor of all Artiodactyla at least 64 million years ago. The amino acid sequence of this gene, termed ARTgagpol, contains recognizable matrix, capsid, nucleocapsid, and reverse transcriptase domains in ruminants, with an RNase H domain in camels and pigs. Phylogenetic analysis and structural prediction of its reverse transcriptase and RNase H domains groups ARTgagpol with gammaretroviruses. Transcriptomic analysis shows ARTgagpol expression in multiple tissues suggestive of a co-opted host function. These findings identify the oldest and largest ERV-derived gagpol gene with an intact ORF in mammals, an intriguing milestone in the co-evolution of mammals and retroviruses. IMPORTANCE Retroviruses are unique among viruses that infect animals as they integrate their reverse-transcribed double-stranded DNA into host chromosomes. When this happens in a germline cell, such as sperm, egg, or their precursors, the integrated retroviral copies can be passed on to the next generation as endogenous retroviruses (ERVs). On rare occasions, the genes of these ERVs can be domesticated by the host. In this study we used computational similarity searches to identify an ancient ERV with an intact viral gagpol gene in the genomes of camels that is also found in the same genomic location in other even-toed ungulates suggesting that it is at least 64 million years old. Broad tissue expression and predicted preservation of the reverse transcriptase fold of this protein suggest that it may be domesticated for a host function. This is the oldest known intact gagpol gene of an ancient retrovirus in mammals.

Genetic diversity and relatedness in captive collared peccaries Dicotyles tajacu (Linnaeus, 1758) (Cetartiodactyla: Tayassuidae) estimated by microsatellite genotyping using high-throughput sequencing: Implications for their conservation and reintroduction.

The global population of Dicotyles tajacu (Linnaeus, 1758) (Cetartiodactyla: Tayassuidae), commonly known as the collared peccary and distributed in the Neotropics, is currently in decline due to anthropogenic pressures. In this study, five microsatellite loci were used to genetically characterize a group of 20 captive-born collared peccaries intended for reintroduction. This study aimed to evaluate the genetic diversity and relatedness of captive individuals using microsatellite markers. The genetic data generated were used to evaluate the viability of the reintroduction and to propose measures for the management and conservation of this species. In this study, we found relatively high genetic diversity indices, indicating that the group was genetically diverse. Inbreeding coefficients with negative values were observed, indicating an excess of alleles in heterozygosis and an absence of inbreeding. One locus showed deviation from Hardy-Weinberg equilibrium, which may have been caused by the mixing of individuals from different origins. Relatedness analysis indicated that some individuals were highly related, with coefficients indicating they may be first-degree relatives. Our findings indicate that the studied group has enough genetic diversity to be released into nature, but the high individual relatedness found would require the adoption of strategies after the release of animals in the wild to ensure their persistence.

Birth-and-death evolution of ribonuclease 9 genes in Cetartiodactyla.

RNase9 plays a reproductive function and has been recognized as an important member of the ribonuclease (RNase) A superfamily, a gene family that is widely used as a model for molecular evolutionary studies. Here, we identified 178 RNase9 genes from 95 Cetartiodactyla species that represent all four lineages and 21 families of this clade. Unexpectedly, RNase9 experienced an evolutionary scenario of "birth and death" in Ruminantia, and expression analyses showed that duplicated RNase9A and RNase9B genes are expressed in reproductive tissues (epididymis, vas deferens or prostate). This expression pattern combined with the estimate that these genes duplicated during the middle Eocene, a time when Ruminantia become a successful lineage, suggests that the RNase9 gene duplication might have been advantageous for promoting sperm motility and male fertility as an adaptation to climate seasonality changes of this period. In contrast, all RNase9 genes were lost in the Cetacean lineage, which might be associated with their high levels of prostatic lesions and lower reproductive rates as adaptations to a fully aquatic environment and a balance to the demands of ocean resources. This study reveals a complex and intriguing evolutionary history and functional divergence for RNase9 in Cetartiodactyla, providing new insights into the evolution of the RNaseA superfamily and molecular mechanisms for organismal adaptations to the environment.

Genomic and anatomical comparisons of skin support independent adaptation to life in water by cetaceans and hippos.

The macroevolutionary transition from terra firma to obligatory inhabitance of the marine hydrosphere has occurred twice in the history of Mammalia: Cetacea and Sirenia. In the case of Cetacea (whales, dolphins, and porpoises), molecular phylogenies provide unambiguous evidence that fully aquatic cetaceans and semiaquatic hippopotamids (hippos) are each other's closest living relatives. Ancestral reconstructions suggest that some adaptations to the aquatic realm evolved in the common ancestor of Cetancodonta (Cetacea + Hippopotamidae). An alternative hypothesis is that these adaptations evolved independently in cetaceans and hippos. Here, we focus on the integumentary system and evaluate these hypotheses by integrating new histological data for cetaceans and hippos, the first genome-scale data for pygmy hippopotamus, and comprehensive genomic screens and molecular evolutionary analyses for protein-coding genes that have been inactivated in hippos and cetaceans. We identified eight skin-related genes that are inactivated in both cetaceans and hippos, including genes that are related to sebaceous glands, hair follicles, and epidermal differentiation. However, none of these genes exhibit inactivating mutations that are shared by cetaceans and hippos. Mean dates for the inactivation of skin genes in these two clades serve as proxies for phenotypic changes and suggest that hair reduction/loss, the loss of sebaceous glands, and changes to the keratinization program occurred ∼16 Ma earlier in cetaceans (∼46.5 Ma) than in hippos (∼30.5 Ma). These results, together with histological differences in the integument and prior analyses of oxygen isotopes from stem hippopotamids ("anthracotheres"), support the hypothesis that aquatic skin adaptations evolved independently in hippos and cetaceans.

The steroid metabolome of pregnancy, insights into the maintenance of pregnancy and evolution of reproductive traits.

Modes of mammalian reproduction are diverse and not always conserved among related species. Progesterone is universally required to supports pregnancy but sites of synthesis and metabolic pathways vary widely. The steroid metabolome of mid-to late gestation was characterized, focusing on 5α-reduced pregnanes in species representing the Perissodactyla, Cetartiodactyla and Carnivora using mass spectrometry. Metabolomes and steroidogenic enzyme ortholog sequences were used in heirarchial analyses. Steroid metabolite profiles were similar within orders, whales within cetartiodactyls for instance, but with notable exceptions such as rhinoceros clustering with goats, and tapirs with pigs. Steroidogenic enzyme sequence clustering reflected expected evolutionary relationships but once again with exceptions. Human sequences (expected outgroups) clustered with perissodactyl CYP11A1, CYP17A1 and SRD5A1 gene orthologues, forming outgroups only for HSD17B1 and SRD5A2. Spotted hyena CYP19A1 clustered within the Perissodactyla, between rhinoceros and equid orthologues, whereas CYP17A1 clustered within the Carnivora. This variability highlights the random adoption of divergent physiological strategies as pregnancy evolved among genetically similar species.

Comparative chromosome painting in Chacoan peccary, Catagonus wagneri.

Chacoan peccary (Catagonus wagneri, 2n=20) is the most endangered of three extant species of Tayassuidae. Its karyotype has been studied only by differential chromosome staining methods so far. To establish a comparative cytogenetic map of the peccary, we used cross-species hybridization with porcine (Sus scrofa, 2n=38) painting probes. Painting revealed 30 evolutionary conserved autosomal segments between pig and peccary. The q-arm of the submetacentric chromosome X is homologous to the porcine X chromosome, while the p-arm is composed of heterochromatin. Nucleolar organizer regions were detected on chromosomes 8 and 9 which are homologous to pig chromosomes 8 and 4/18, respectively. Fusions of chromosomes homologous to pig chromosomes 4/7 and 4/18 and fission of chromosome 7 are synapomorphic characters shared by Catagonus wagneri and Tayassu pecari but not by Pecari tajacu. Our results confirmed a high rate of karyotype evolution in Tayassuidae and a closer relationship of Catagonus wagneri with Tayassu pecari than with Pecari tajacu.

Reassessment of the enigmatic ruminant Miocene genus Amphimoschus Bourgeois, 1873 (Mammalia, Artiodactyla, Pecora).

Amphimoschus is an extinct Eurasian ruminant genus, mostly recorded in Europe, without a close living relative and, hence, an unknown systematic position. This genus is known from around 50 localities from the late early to the middle Miocene. Two species were described during 180 years, but since their first description during the late 19th century and early 20th century, hardly any detailed taxonomic work has been done on the genus. Over the years, extensive collecting and excavating activities have enriched collections with more and more complete material of this still rare and enigmatic animal. Most interestingly, a number of skull remains have been unearthed and are promising in terms of providing phylogenetic information. In the present paper, we describe cranial material, the bony labyrinth, the dentition through 780 teeth and five skulls from different ontogenetic stages. We cannot find a clear morphometric distinction between the supposedly smaller and older species Amphimoschus artenensis and the supposedly younger and larger species A. ponteleviensis. Accordingly, we have no reason to retain the two species and propose, following the principle of priority (ICZN chapter 6 article 23), that only A. ponteleviensis Bourgeois, 1873 is valid. Our studies on the ontogenetic variation of Amphimoschus does reveal that the sagittal crest may increase in size and a supraorbital ridge may appear with age. Despite the abundant material, the family affiliation is still uncertain.

Sequencing and characterization of the complete mitochondrial genome of Mishmi takin (Budorcas taxicolor taxicolor) and comparison with the other Caprinae species.

Mishmi takin (Budorcas taxicolor taxicolor) is an endangered and least studies species. We sequenced and annotated a 16,584 bp long complete mitogenome of Mishmi takin (B. t. taxicolor). It was encoded by 37 genes: 13 Protein-Coding Genes (PCGs), two ribosomal RNA genes, 22 transfer RNA, and non-coding control region. All tRNA genes show a typical secondary cloverleaf structure, except tRNAser(AGY) in which dihydrouridine arm did not form a stable structure. We observed 82 bp nucleotide deletions in the control region of Mishmi takin. The overall analysis of 13 PCGs, two rRNA and 22 tRNA of Mishmi takin indicates its close relationship with the B. t. tibetana than the B. t. bedfordi. The overall genetic similarity indicated that all takin species is closely related to the P. nayaur than the C. ibex and Ovibos moschatus. This study provides a useful resource for further studies on conservation and population genetics of Budorcas.

The complete mitochondrial genome of Budorcas taxicolor tibetana (Artiodactyla: Bovidae) and comparison with other Caprinae species: Insight into the phylogeny of the genus Budorcas.

Takin (Budorcas taxicolor) is one of the most endangered species. However, the taxonomy of takin is still in dispute. Here, the complete mitochondrial genome of Budorcas taxicolor tibetana was reported for the first time, which featured a typical circular molecule of 16,665 bp in length, including 13 PCGs, 22 tRNAs, two rRNAs, and one control region. A + T content was higher than G + C content. All of the genes were encoded on the heavy strand, except for eight tRNAs and ND6 gene. The OL region was 49 bp in length and highly conserved in the synthesis and stem-loop regions, and all of the substitutions and indels were found only in the loop structure. Three types of tandem repeat units, six pairs of hairpin loop structure (TACAT, ATGTA) and six CSBs were predicted in the control region. Our results clearly revealed the systematic status of Budorcas species, and the phylogenetic analyses indicated that Budorcas was closer to Capra and Pseudois, rather than to Ovis, which should be merged into the subfamily of Caprinae.

Signature of positive selection in mitochondrial DNA in Cetartiodactyla.

Acceleration of the amino acid substitution rate is a good indicator of positive selection in adaptive evolutionary changes of functional genes. Genomic information about mammals has become readily available in recent years, as many researchers have attempted to clarify the adaptive evolution of mammals by examining evolutionary rate change based on multiple loci. The order Cetartiodactyla (Artiodactyla and Cetacea) is one of the most diverse orders of mammals. Species in this order are found throughout all continents and seas, except Antarctica, and they exhibit wide variation in morphology and habitat. Here, we focused on the metabolism-related genes of mitochondrial DNA (mtDNA) in species of the order Cetartiodactyla using 191 mtDNA sequences available in databases. Based on comparisons of the dN/dS ratio (ω) in 12 protein-coding genes, ATP8 was shown to have a higher ω value (ω = 0.247) throughout Cetartiodactyla than the other 11 genes (ω < 0.05). In a branch-site analysis of ATP8 sequences, a markedly higher ω value of 0.801 was observed in the ancestral lineage of the clade of Cetacea, which is indicative of adaptive evolution. Through efforts to detect positively selected amino acids, codon positions 52 and 54 of ATP8 were shown to have experienced positive selective pressure during the course of evolution; multiple substitutions have occurred at these sites throughout the cetacean lineage. At position 52, glutamic acid was replaced with asparagine, and, at position 54, lysine was replaced with non-charged amino acids. These sites are conserved in most Artiodactyla. These results imply that the ancestor of cetaceans underwent accelerated amino acid changes in ATP8 and replacements at codons 52 and 54, which adjusted metabolism to adapt to the marine environment.

Characterization of the complete mitogenome of Indian Mouse Deer, Moschiola indica (Artiodactyla: Tragulidae) and its evolutionary significance.

The mitochondrial genome of Indian mouse deer (Moschiola indica) was sequenced, assembled and characterized for the first time using 22 pairs of polymerase chain reaction (PCR) primers. The mitogenome of M. indica which is 16,444 bp in size was found very similar to most vertebrates in organisation that harbours 13 protein-coding genes, 22 transfer RNA, 2 ribosomal RNA and 1A + T-rich region. Its comparison with over 52 mitogenomes of the order Artiodactyla, showed a conserved nature of gene organisation, codon usage, gene orientation and evolutionary rates of proteins except that M. indica possesses an extra copy of trnF. The complete mitogenome and protein-coding genes of M. indica were found to be highly A + T biased. Rate of protein evolution was highest in atp8 and lowest in cox3. Further, a higher purifying selection pressure was found to be acting on family Tragulidae compared to Bovidae and Cervidae. The phylogenetic analysis of M. indica placed the Tragulidae as sister-group of all other ruminants, similar to previous analyses. Moschiola forms the sister-group to the other two tragulid genera Tragulus (from Asia) and Hyemoschus (from Africa), which is unexpected as usually the Asian species are thought to form a monophyletic group.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing.

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

Ancient DNA analysis of the extinct North American flat-headed peccary (Platygonus compressus).

The geographical range of extant peccaries extends from the southwestern United States through Central America and into northern Argentina. However, from the Miocene until the Pleistocene now-extinct peccary species inhabited the entirety of North America. Relationships among the living and extinct species have long been contentious. Similarly, how and when peccaries moved from North to South America is unclear. The North American flat-headed peccary (Platygonus compressus) became extinct at the end of the Pleistocene and is one of the most abundant subfossil taxa found in North America, yet despite this extensive fossil record its phylogenetic position has not been resolved. This study is the first to present DNA data from the flat-headed peccary and full mitochondrial genome sequences of all the extant peccary species. We performed a molecular phylogenetic analysis to determine the relationships among ancient and extant peccary species. Our results suggested that the flat-headed peccary is sister-taxon to a clade comprising the extant peccary species. Divergence date estimates from our molecular dating analyses suggest that if extant peccary diversification occurred in South America then their common ancestor must have dispersed from North America to South America well before the establishment of the Isthmus of Panama. We also investigated the genetic diversity of the flat-headed peccary by performing a preliminary population study on specimens from Sheriden Cave, Ohio. Flat-headed peccaries from Sheriden Cave appear to be genetically diverse and show no signature of population decline prior to extinction. Including additional extinct Pleistocene peccary species in future phylogenetic analyses will further clarify peccary evolution.

Development of deft amplification refractory mutation sequencing system (ARMSS) for discriminating Pilos antler based on a short cytochrome b (Cytb) gene.

Pilos antler (Lu-Rong in Chinese) is a famous traditional medicine in China. Many adulterants have been discovered in Chinese markets in recent years. However, few DNA-based methods are effective for discrimination of this DNA-degraded animal medicine. Here, novel and deft amplification refractory mutation sequencing system (ARMSS), integrating the advantages of the amplification refractory mutation system (ARMS) and the short DNA barcode, was first developed to discriminate Pilos antler from its adulterants. We aimed to provide a new sight and inspiration for deft detection. The results showed that developed ARMS achieved strong specificity and high sensitivity in rapid identification, while the short Cytb gene was of excellent identification power in terms of accurate identification, which suggested that ARMSS successfully integrated the advantages of the ARMS and short DNA barcode, and that it was useful for deft detection. Our study determined that the deft ARMSS could be the well candidate for discrimination of Pilos antler, as well as be a valuable tool for deft identification of Chinese medicine.

Inbreeding and Offspring Sex Ratio in the Pygmy Hippopotamus (Cheoropsis liberiensis) Population Kept in Zoological Gardens.

The aim of this study was to estimate the inbreeding level and its impact on offspring sex ratio in the pygmy hippopotamus population kept in zoological gardens. Records of pygmy hippopotamus born between 1873-2013 were extracted from the international studbook. Totally, 1357 individuals originating from 148 breeding units were included (individuals with unknown sex were omitted). The offspring sex ratio is defined as the number of sons to the total number of progeny of each dam and sire. Spearman's rank correlation was employed to examine the relationships between the inbreeding level and offspring sex ratio. Inbreeding coefficients and individual increase in inbreeding coefficients (included as a linear co-variable) were examined as well as the geographic region and birth period using general linear models. The average inbreeding coefficient was 5.39%. The following sex proportion was observed for the inbred population: 57% and 43% for females and males, respectively. A significant relationship between inbreeding level of parents and their offspring sex ratio were estimated for European zoological gardens, whereas in others geographic regions the dependencies were insignificant.

Placental PAGs: gene origins, expression patterns, and use as markers of pregnancy.

Pregnancy-associated glycoproteins (PAGs) are abundantly expressed products of the placenta of species within the Cetartiodactyla order (even-toed ungulates). They are restricted to this order and they are particularly numerous in the Bovidae. The PAGs exhibit a range of temporal and spatial expression patterns by the placental trophoblasts and probably represent a group of related proteins that perform a range of distinct functions in the epitheliochorial and synepitheliochorial placental forms. This review presents an overview of the origins of the PAGs, a summary of PAG expression patterns, and their use as markers of pregnancy status. Speculations about their putative role(s) in pregnancy are also presented.

Hippos stem from the longest sequence of terrestrial cetartiodactyl evolution in Africa.

According to molecular data, hippopotamuses and cetaceans form a clade excluding other extant cetartiodactyls. Despite a wealth of spectacular specimens documenting cetacean evolution, this relationship remains poorly substantiated by the fossil record. Indeed, the evolutionary path leading from the hippo-cetacean ancestor to Hippopotamidae is plagued by missing fossil data and phylogenetic uncertainties. Only an origination within the extinct anthracotheres is compatible with molecular results, substantial filling of phyletic gaps and recent discoveries of early Miocene hippopotamids. Yet, the anthracothere stock that gave rise to Hippopotamidae has not been identified. Consequently, recent phylogenetic accounts do not properly integrate the anthracotheriid hypothesis, and relate Hippopotamidae to a stretched ghost lineage and/or close to Suina. Here we describe a new anthracothere from Lokone (Kenya) that unambiguously roots the Hippopotamidae into a well-identified group of bothriodontines, the first large mammals to invade Africa. The hippos are deeply anchored into the African Paleogene.

Reliable discrimination of 10 ungulate species using high resolution melting analysis of faecal DNA.

Identifying species occupying an area is essential for many ecological and conservation studies. Faecal DNA is a potentially powerful method for identifying cryptic mammalian species. In New Zealand, 10 species of ungulate (Order: Artiodactyla) have established wild populations and are managed as pests because of their impacts on native ecosystems. However, identifying the ungulate species present within a management area based on pellet morphology is unreliable. We present a method that enables reliable identification of 10 ungulate species (red deer, sika deer, rusa deer, fallow deer, sambar deer, white-tailed deer, Himalayan tahr, Alpine chamois, feral sheep, and feral goat) from swabs of faecal pellets. A high resolution melting (HRM) assay, targeting a fragment of the 12S rRNA gene, was developed. Species-specific primers were designed and combined in a multiplex PCR resulting in fragments of different length and therefore different melting behaviour for each species. The method was developed using tissue from each of the 10 species, and was validated in blind trials. Our protocol enabled species to be determined for 94% of faecal pellet swabs collected during routine monitoring by the New Zealand Department of Conservation. Our HRM method enables high-throughput and cost-effective species identification from low DNA template samples, and could readily be adapted to discriminate other mammalian species from faecal DNA.

Genomic organization and differential signature of positive selection in the alpha and beta globin gene clusters in two cetacean species.

The hemoglobin of jawed vertebrates is a heterotetramer protein that contains two α- and two ß-chains, which are encoded by members of α- and ß-globin gene families. Given the hemoglobin role in mediating an adaptive response to chronic hypoxia, it is likely that this molecule may have experienced a selective pressure during the evolution of cetaceans, which have to deal with hypoxia tolerance during prolonged diving. This selective pressure could have generated a complex history of gene turnover in these clusters and/or changes in protein structure themselves. Accordingly, we aimed to characterize the genomic organization of α- and ß-globin gene clusters in two cetacean species and to detect a possible role of positive selection on them using a phylogenetic framework. Maximum likelihood and Bayesian phylogeny reconstructions revealed that both cetacean species had retained a similar complement of putatively functional genes. For the α-globin gene cluster, the killer whale presents a complement of genes composed of HBZ, HBK, and two functional copies of HBA and HBQ genes, whereas the dolphin possesses HBZ, HBK, HBA and HBQ genes, and one HBA pseudogene. For the ß-globin gene cluster, both species retained a complement of four genes, two early expressed genes-HBE and HBH-and two adult expressed genes-HBD and HBB. Our natural selection analysis detected two positively selected sites in the HBB gene (56 and 62) and four in HBA (15, 21, 49, 120). Interestingly, only the genes that are expressed during the adulthood showed the signature of positive selection.

Heterogeneous tempo and mode of conserved noncoding sequence evolution among four mammalian orders.

Conserved noncoding sequences (CNSs) of vertebrates are considered to be closely linked with protein-coding gene regulatory functions. We examined the abundance and genomic distribution of CNSs in four mammalian orders: primates, rodents, carnivores, and cetartiodactyls. We defined the two thresholds for CNS using conservation level of coding genes; using all the three coding positions and using only first and second codon positions. The abundance of CNSs varied among lineages, with primates and rodents having highest and lowest number of CNSs, respectively, whereas carnivores and cetartiodactyls had intermediate values. These CNSs cover 1.3-5.5% of the mammalian genomes and have signatures of selective constraints that are stronger in more ancestral than the recent ones. Evolution of new CNSs as well as retention of ancestral CNSs contribute to the differences in abundance. The genomic distribution of CNSs is dynamic with higher proportions of rodent and primate CNSs located in the introns compared with carnivores and cetartiodactyls. In fact, 19% of orthologous single-copy CNSs between human and dog are located in different genomic regions. If CNSs can be considered as candidates of gene expression regulatory sequences, heterogeneity of CNSs among the four mammalian orders may have played an important role in creating the order-specific phenotypes. Fewer CNSs in rodents suggest that rodent diversity is related to lower regulatory conservation. With CNSs shown to cluster around genes involved in nervous systems and the higher number of primate CNSs, our result suggests that CNSs may be involved in the higher complexity of the primate nervous system.